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Abstract 

The mass transfer of interstitial impurities in a crystalline lattice under the influence 
of the fast-moving deformation disturbance of the type of a shock wave is considered. 
The velocity of the movement of the disturbance is supposed to be compared with 
the characteristic velocity of the relaxation of the diffusion flux to its local equilib- 
rium value determined by the Pick's law. The similar situation occurs in a number 
of experiments on the exposure of a solid to dynamical external loads giving rise to 
such fast hydrodynamical processes in a sample that the local equilibrium assump- 
tion, normally assumed for the macroscopic description of transport processes, is no 
longer valid. Considering the diffusion flux among the set of independent variables 
we have derived a set of coupled hydrodynamic equations describing nonequilibrium 
behavior of a solid in the absence of local equilibrium in the system. Within the 
scope of the proposed model it has been shown that in comparison with the local 
equilibrium system an enhanced mass transfer occurs under local nonequilibrium 
conditions. 
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NOMENCLATURE 



c mass concentration of interstitialcs 
E energy density 

J diffusion flux density of intersti- 
tials 

p pressure 

Q energy flux density 

q heat flux density 

S entropy density of the medium 

T temperature 

t time 

u displacement vector 
Uij infinitesimal strain tensor 
V velocity of the deformation distur- 
bance 

Vd speed of the propagation of con- 
centration perturbations 



V mass velocity 
Xi coordinate 

Greek symbols 

5 width of the transition region 

// chemical potential of the medium 

V chemical potential of the intersti- 
tials per unit of volume 

^ dimensionless coordinate 
liik the tensor momentum flux den- 
sity 

TTjj viscous pressure tensor 

p total mass density 

(Tij elastic stress tensor 

T relaxation time of the diffusion flux 

if =V/Vd constant value 
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1 Introduction 



In recent experiments [1,2] on the irradiation of the thin- film metal samples by 
high power ion beams (HPIB) thin films of aluminium and silver on the nio- 
bium and copper substrates were irradiated by carbon ions. The considerable 
penetration depth of the film atoms into the substrate, exceeding the length 
of the projective path of carbon atoms, has been revealed in the irradiated 
samples. The estimation of the effective diffusion coefficients has given the 
anomalously high magnitudes [1-3] D ~ 10~^ — 10~^ sm^s~^. The effect of the 
intensive mass transfer by HPIB irradiation has been registered earlier in work 
[4] . A similar effect has also been observed at electron and laser irradiation of 
solids [5-7] and at shock loading of the samples [8,9]. The complete explana- 
tion of this phenomenon within the scope of the known diffusion mechanisms 
has not yet been found. 

One of possible causes responsible for the accelerated mass transport is the 
shock wave, generated by irradiation, which is able to induce the enhanced 
migration of atoms on large depths [2]. In investigating this mechanism of 
diffusion within the scope of the hydrodynamic description usually one starts 
from the known diffusion equation in which, according to Pick's law, the diffu- 
sion flux is deflned by gradients of the concentration of a diffusing substance, 
a pressure and a temperature. Such approach is well justified if the distur- 
bances occurring in a system change slowly enough in comparison with the 
characteristic rate at which the given system approaches the local equilibrium. 
Meanwhile the interaction of particle beams with a solid induces in the latter 
rather fast hydrodynamic processes with the characteristic times of order or 
less than lO^^s and the shock velocity can reach several kilometers per secund. 
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Under these conditions one should expect a significant deviation of the sys- 
tem state from local equilibrium and the violation of Pick's law and take into 
account an influence of the relaxation of the system on the process of mass 
transport. 

In this paper we want to study within a simple model basic features the im- 
purity transfer in the relaxating solid medium. More specifically the purpose 
of the work is to research the interstitial diffusion in the field of the fast mov- 
ing deformation disturbance of the type of the shock wave. The velocity of 
the movement of this disturbance is supposed to be compared with the ve- 
locity of the propagation of concentration disturbances in the medium. The 
latter condition just defines the condition of absence of local equilibrium in 
the system [11]. The task is also interest in connection with shock-compression 
experiments with metals containing considerable atomic concentrations of in- 
terstitial hydrogen [12,13]. 

The description of the fast hydrodynamical processes requires an extended 
number of independent variables, defining the nonequilibrium state of a medium, 
in comparison with conventional hydrodynamics. During the course of short 
enough the characteristic time intervals, in particular, the diffusion fiux may 
not have time to relax to its local-equilibrium value, defined by Pick's law 
and it should be considered as independent variable obeying some relaxation 
equation. The dissipative fiuxes (diffusion fiux, viscous pressure tensor, heat 
fiux) included in the local conservative laws of mass, momentum and energy 
are frequently exploited as additional independent variables [10,11]. 
In section 2 we determine the total set of the independent variables suitable 
for the description of the fast hydrodynamical processes in the crystalline solid 
with interstitials. Por these variables the derivation of the set of the coupled 
hydrodynamical equations, including a relaxation equation for the diffusion 
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flux is given. Section 3 is devoted to the investigation of the mass transport in 
the fleld of fast-moving deformation disturbance of the crystal in the absence 
of local equilibrium in the system. Resume and concluding remarks are given 
in section 4. 



2 Derivation of hydrodynamic equations 



In this section we want to derive the set of equations relevant for the descrip- 
tion of fast hydrodynamical processes in the crystalline solid with interstitial 
impurities. The complete set of the hydrodynamical equations has to contain 
the local conservation laws of mass , momentum and the equation for the 
entropy production 

|^+divpv = 0, (1) 

— + div (,5v + - — J = - . (3) 

Here p is the total mass density of the medium (lattice plus interstitials) , v 
is the mass velocity, S is the entropy density of the medium, Iljfc is the ten- 
sor momentum flux density, q is the heat flux density, J is the diffusion flux 
density of interstitials, v and T are the chemical potential of the interstitials 
per unit of volume and the absolute temperature, respectively, and R{R> 0) 
is the dissipative function of the medium; summation over repeated Latin in- 
dices is implied. 

Further we suppose that a concentration of vacancies is negligible or their 
mobility is signiflcantly less than the mobility of interstitials so that the mass 
transport occurs mainly by the interstitial motion. Introducting the mass con- 
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centration of the interstitial particles c we write down the continuity equation 
for the interstitial density cp as follows 



^ + div(cpv + J) = 0. (4) 



It should be noted that the diffusion flux J, appearing in equations (3) and (4), 
is determined in a rest frame (moving with velocity v relative to the labora- 
tory frame). Meanwhile for a crystal it is convenient to determine the diffusion 
flux relative to a lattice. It is easy to see that both fluxes are related by 
J = (1 — c)Jj^. As the mass concentration c ~ mp/m^, where rUp and are 
masses of an impurity particle and a host atom, then for light impurities (type 
of hydrogen), occurring into a matrix of heavy atoms, c << 1 even for the 
finite atomic concentration. Because of this J ~ J^. 

The variables p, v, c, (or T) define conventional set of the independent hy- 
drodynamical variables for a "fluid-hke" medium. A solid with point defects 
has the extra independent hydrodynamical variable. This is a displacement 
vector u defining displacement of lattice sites [14] . In addition, according to 
above said, in the case of fast processes the dissipative fluxes - diffusion flux 
J, heat flux q and viscous pressure tensor TTjfe must be included among the set 
of independent variables. For simplicity we shall assume that the relaxation 
times of fluxes q and TTjfc are considerably less than the relaxation time of dif- 
fusion flux and during characteristic time of the process under consideration 
the fluxes q and Trjfe have time to relax to their local equilibrium values [11]. 
In this case diffusion flux J is the only new independent variable. 
Thus we assume that quantities p, v, c, u, J completely deflne the nonequi- 
hbrium state of the system. For such set variables the set of the hydrody- 
namical equations (l)-(4) has to be supplemented by the equations for the 
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displacement vector u and the diffusion flux J. To derive these equations and 
define the still unknown Uik, q and R, appearing in equations (l)-(4), we use 
the method originally applied in the theory of the superfluids for the derivation 
of the equations of the two-fluid hydrodynamics [15]. The main idea, in the 
given case, is as follows. We shall find these quantities so that the conservation 
energy law 

OE 

— + divQ = 0, (5) 

would follow from equations (l)-(4) and the equations for u, J (an upper dot 

denotes the time derivative). In the equation (5) E and Q are the energy 

density and the energy flux density of the medium, respectively. 

The energy E is related by the Gallilean transformation to its value Eq in 

the frame, where a given element of the volume of the medium rests, by the 

relationship 

E^Eo + ^. (6) 

Let us now write down the differential of Eq, considered as a function of 
S, p, c, J and the infinitesimal strain tensor Uij in the form 

dEo — TdS + jidp + vdc + dijduij + w • d J , (7) 

here p is the chemical potential of the medium, w = dEo/dJ is the conjugate 
variables with J, and (7^ is the symmetric, elastic stress tensor; here and 
further the tilde is used to denote the traceless part of a tensor. In the equation 
(7) it has been taken into account that the variation duu of sum of diagonal 
components of the strain tensor Uij is determined by the variation of the 
density dp. 
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Differentiating the equation (6) with respect to time and using equation (7), 
one obtains 



dE d /pf \ ^dS dp dc dui ■ , , 



Employing equations (l)-(4) we can write down the time derivatives as 



di^ ll 2 dt ' dt 

= -div(^v) - Vi — {Uik - pViVk) , (9) 

r^ = -rdiv(5v + |--^j)+i? 

= -div(r5v + q - - j) + V • V{TS) -Tv-VS 

_J.v(^)+« + C.^+TJ.vQ, (10) 

i/c= -i/v Vc-div(^j) + J- V(^) , (11) 
HP = -div(//pv) + V • V(/ip) - /iv • Vp , (12) 



a. 



diii d . . .da, 



dijiii) - Ui— — . (13) 



dxj dxj ^" dxj 



Inserting the equations (9)-(13) into equation(8) and taking into account that 



VEo = TVS + hVp + i/Vc + dijV [—^ ) + wyJi , 
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after long, though simple, transformations, one has 



E = -div{^v + TSw + /xpv + q - li • a} + ^fe^ii^(^) 
d 

(Iljfc - pViVk + CTik) + V • V(--E'o + TS + (ip) 



dxk 



Ui 



T \pT^ ' " dxk 

+Wi{ji + v -VJi) (14) 



Finally equation (14) can be rewritten as 
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div{^v + TS\ + /xpv + q- u- a + v-7r} 



Vi avi / 1/ N dui aaik aji 

-R + ci—+ vr..— + TJ . V(-) + - + , (15) 

where the viscous pressure tensor iTik is defined from the equahty 

Hife = pViVk + p5ik - + -Kik (16) 

and the notations have been introduced 

p=-Eo + TS + lip, (17) 

and for the material derivative d/dt — d/dt + v - V . 

In the derivation equation (16) we have neglected the term of the second order 
of infinitesimal in strains djkdui/ dxj in comparison with the linear term dik- 
Equation (16) defines the momentum fiux density, where p can be interpreted 
as "pressure" in the medium. Comparing further equation (15) with the energy 
conservation law (5) we define the energy fiux density Q and the dissipative 
function R as 

2 

Q= (^ + So+p)v + q-u-CT + v-7r (18) 



and 



Since the energy Eq is even under time reversal, the dependence between Eq 
and Ji in the simplest approximation has the form 

^0 = ^ WijJiJj + . . . , Wij^ Wji , 
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where Wij is assumed to be constants and dots denote the terms of the highest 
order of infinitesimal in fluxes. In this case 

OE 

= ^ = (20) 
oJi 

and the dissipative function R can be rewritten as 

VT dvi dui daik 

T oxk dt oxk 

where 

dxi^pT^ dt 

In the framework of the hnear theory the positive definiteness of R leads to 
the linear relationships relating the dissipative fluxes to the thermodynamic 
forces. Taking into account Onsager's reciprocity relations for the transport 
coefficients and time-reversal property of dissipative effects [16], we can write 
these relationships in the form 



T^ik = -Viklm-K , [^^) 

OXffi 

„ _ <k dT , , dakj , . 

- oiik^k - Pik^— > y^'^) 



dxk " ^^'^ dxj ' 

T <k 9T , dakj 

dt~ T^dx, ^''dxj' ^^^^ 

where the forth-rank tensor rj is related to the effect of the viscosity and 
a[j^, Plj^, . . . are tensor transport coefficients. 

Let us eliminate the thermodynamic force from equations (23) and (25). 
Redenoting the transport coefficients and introducing the explicit definition of 
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Xfe into equation (24) one rewrites the set (23)-(25) for the case of the cubic 
lattice as 




VT 

q — aJ — — — /?V • a , 



(26) 



(27) 



-XV -a, 



(28) 



where a, /3, 7 . . . are the scalar transport coefficients, k is connected with the 
heat conductivity , r is the relaxation time of the diffusion flux. The coefficients 
K, A determine the relation between the dissipative effects and shear stresses. 
The remaining coefficients are due to the processes of diffusion, thermal dif- 
fusion and barodiffusion (owing to the pressure dependence of the interstitial 
chemical potential 1/ in equation (28)). 

Thus the complete set of the hydrodynamic equations of the system under con- 
sideration contains now equations (1)- (4) and equations (27), (28), where the 
fluxes 7rjfe,q and the dissipative function R are given by equations (22), (26) 
and (21). As it seen from (28) the equation for diffusion flux is the equation 
of a relaxation type. At r = 0, i.e. for slow processes, it determines directly 
the diffusion flux in terms of gradients of the basic hydrodynamic variables. 

3 Mass transfer in field of fast-moving deformation disturbance 

In this section we consider the diffusion of interstitial impurities in the field 
of the fast-moving deformation disturbance of the type of a chock wave. We 

assume that the stepped disturbance (kink), for which the displacement field 
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(along the x-axis) has the form [17,18] 

x-Vt 



X, t) = Uo(l- tanh — - — ) , (29) 



where uq is the amphtude of the disturbance and 5 is the width of the tran- 
sition region, propagates along the x-direction with the constant velocity V 
compared with the velocity of the propagation of concentration disturbances 
in the medium. 

At the given motion of the lattice (29) the diffusion process is largely governed 
by the set of two equations: the local conservation law of mass (4) and the 
relaxation equation for the diffusion flux (28). Let us now transform the right 
side of equation (28) to the alternative form. Combining equations (7) and 
(17) and taking into account equation (20), one obtains 

W 

pdjjt, — —SdT + dp + ud c + aijduij + —dJ , (30) 



where W — ^Wu. It follows from equation (30) that one can consider the chem- 
ical potential to be a function of the independent variables T,p,c,Uij, J^. 
Omitting essentially nonlinear terms, note that in the case of the cubic sym- 
metry a scalar function can depend on the tensor Uij by convolution uu — 
only. Consequently one can consider u to he a function of T, p, c only. 
The available experimental data indicate that when a shock wave propagates in 
a a solid the effect of thermodiffusion is neglected in comparison with borod- 
iffusion [2,19]. Here we also assume that the main factor determining mass 
transfer in our system is the effect of barodiffusion. Then omitting the gra- 
dient of a temperature one rewrites equation (28) in the more simple form 
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as 



+ J = - 7T#(^) Vc - 7^(^) Vp - AV • a , (31) 



dt dc^pT' dp^pT 



It follows from equation (16) that = —pSik + d'ik, with era = —3p is the 
elastic stress tensor. In the one-dimensional case one neglects the influence 
of shear stress on a diffusion mobility Q. In addition, we assume that Hook's 
law p = —kU, relating the pressure p with the deformation U = du/dx of 
a material element, is valid. Redenoting the constants in equation (31) one 
writes down in the linear approximation the coupled set of the one-dimensional 
versions of equations (4) and (31) as 

dc dJ , , 

where G > has the meaning of the coefficient of the barodiffusion and 
D = '-)p^^Td{v / pT) / dc is the diffusion constant. In addition, in obtaining 
of equation (32) and (33) we have neglected the term v ■ V, appearing in the 
material derivative that is small for fast processes in comparison with the time 
derivative d/dt. 

After some algebra, equations (32) and (33) are brought to the form 
dc D d^c _ d'^c _ d^U 

dt ^ V3 dt^ dx^ ~ ^ dxdt • ^ ' 



^ In the one-dimensional case a is proportional to du/dx. Therefore allowance for 
the term AV • a in the right-hand side of Eq.(31) leads simply to redefinition of the 
constant G in Eq. (33). 
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At U — equation (34) (as well as equation(35)) is known as the telegrapher 
equation that, unlike to the classical diffusion equation, gives rise to the finite 
speed = D/t of the propagation of concentration perturbations [10,11]. 
It is convenient to consider equations (34) and (35) in the reference frame 
moving with velocity V together with the front of the disturbance (29). Passing 
to the dimensionless coordinate i — {x — Vt)/TVD, one has 

2 ,.d^c dc d'^lJ 
2 ..d'^J dJ d^U 

where (p = V/Vo, U — GU/cqD, c = (c — Co)/cq, J = J/cqpVd and Cq is the 
equilibrium interstitial concentration. 

Supposing that far from the front of the wave (^ = 0) the medium is not 

disturbed, we take the boundary conditions in the form 

^^±oo, c^O, J^O. (38) 



The solution of the boundary value problem (36) - (38) can be presented as 
follows 



hU{i)-ab J e-<^-^')U{i')di',^<l, 

— oo 

= 1, 



dXJ_ 
'di 



(39) 



hU{0+cihje-<^-^'m{^')d^', (^>1, 



Here a — (pb, b — 1/(1 — ip^) and the solution for the dimensionless fiux J 
follows from (39) as a result of the replacement of b by a. Equation (39) is the 
solution of boundary value problem for arbitrary deformation disturbance U , 
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vanishing at ^ — > ±00. 

The transition to the local-equihbrium theory occurs if in equation (33) one 
puts formally r = 0. Using the same dimensionless variables as in equations 
(36) and (37) one has in this case the following equations 





dCieq 


d^U 






de ' 


d^Jleq 


dJleq 


d^U 


de 


di - 


-^de 



(40) 
(41) 

The solution of this equations has the standard form and we shall not give it 
here. 

For the disturbance (29) U{^) — —Uo/{cosh'^^/6i), where Uq is the known 
constant and di = S/tVd, one rewrites equation (39) in the form 

i 



—00 



cosh ^ 

2 sinh ^ 
Si cosh^ ^ ' 

b 

cosh^ ^ 



— d^,ip<l, 



cosh ^ 



^ = 1, 



e-a<5i(|-«) 

{a6i)b I ,2. 95 >1, 



(42) 



cosh i 



here ^ = i/5i. Equation (42) is the solution of the given problem at arbi- 
trary values of parameters. Here we are interested in values of the parameter 
if — V/Vn close to one for which \a\ » 1, \b\ » 1. At the same time the 
parameter 5i, characterizing the width of the transition region of the distur- 
bance, can take, generally speaking, any values. In what follows we confine 
ourselves to two hmiting cases. 

l.The wide transition region, Si — S/tVe, ~ 1. Figures 1-3 present the curves 
for the impurity concentration c/Uq (as well as for the diffusion flux J/(pUo) 
at different values of the parameter if in the case of the propagating local 
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compression {Uq > 0). This figures also show the difference between the local 
equilibrium and the local nonequilibrium (presence of the relaxation process 
in the medium) theories. In all cases the deformation disturbance produces the 
enhanced impurity concentration immediately before the front of the defor- 
mation wave (the region of > 0), lowering it behind the front. In the region 
of the of the enhanced concentration the diffusion flux points in the direction 
of the disturbance propagation, while in the region of lowered concentration - 
in the opposite direction. 

As figures 2 and 3 show, at </? ~ 1 the concentration perturbations are signif- 
icantly greatly in the local nonequilibrium medium. In this case the diffusion 
flux is nearly twice as large as that for the local equilibrium system. 
2. The narrow transition region, 5i — 5/tVd « 1- In what follows we con- 
sider the special case oi 5i = l/|a|, \a\ » 1. Figures 4-5 demonstrate 
a pronounced difference between two approaches. The ejection of the major 
number of interstitials into the region before the front of the wave is observed 
in the absence of a local equilibrium in the medium. In the latter case the 
impurity concentration in the indicated region is nearly by an order greater 
than the local equilibrium concentration that does not change practically with 
the rising of ip. Figures 4-5 also show the considerable deviation for diffusion 
fluxes demonstrating the acceleration of the mass transfer under local nonequi- 
librium conditions. Note that at > 1 in the region of the front the diffusion 
fluxes referring to the different nonequilibrium states of the medium point in 
complete opposite directions. The difference between the local nonequilibrium 
curves at different magnitudes (p are shown in flgure 6. 

The reason of the enhanced impurity concentration before the front in the 
local-nonequilibrium system is that when the velocity of the deformation wave 
is near to Vd, the concentration disturbances have no time to propagate on 
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large distances from the front. This has the effect of accumulation of intersti- 
tials in the indicated region. Here is a complete analogy with the task about 
the motion of a thermal source in a relaxating medium [20] . 
It should be noted that in the case of the local tension (Uq < 0) the concen- 
tration peaks will form behind the front of the wave < 0) by the moving 
particles from the region ^ > 0. 

4 Summary and conclusion 

In this work we have presented the investigation of mass transfer in the field 
of the fast-moving deformation disturbance of the type of a shock wave, prop- 
agating in the infinite lattice with interstitial impurities. The research is mo- 
tivated by experiments on the exposure of solid to dynamical external loads 
[1-9] 

For the description of fast hydrodynamical processes the extending of the set 
of the conventional variables, characterizing the nonequilibrium state of a solid 
is needed. In the given model the only extra variable is the diffusion fiux, the 
relaxation time of which is assumed to be comparable to the timescale of the 
process under consideration. For the selected variables the set of the coupled 
hydrodynamical equations has been derived. 

Having assumed that the main factor impacting on mass transfer in our case 
is the effect of barodiffusion [2,19], we reduce the diffusion task to the in- 
vestigation of the set of two linear, one-dimensional equations. These are the 
equation for the interstitial concentration (36) and the relaxation equation for 
the diffusion fiux (37). The solution of this set is defined by two dimensionless 
parameters if = V/Vd and = 6/tVd- At the different magnitudes 9? ~ 1 we 
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have considered two special cases 5i ~ 1 and 5i « 1. 

Figures 1-5 show that the occurrence of the relaxation process in the medium 
leads to the significant deviation of the concentration profiles and the diffusion 
fluxes from the local equilibrium diffusion data. The ejection of interstitials 
into the region before the front enables to speak about the effect of the en- 
trainment of impurity particles by the wave. This effect takes place in the local 
equilibrium system also (see [21] as well), however at 5i << 1 to a largest mea- 
sure it is displayed in the local nonequilibrium conditions (figures 4 and 5). 
In the case of the local tension (C/q < 0; the mirror reflection of the plots 
1-6 relative to axis ^) the concentration peaks will form beyond the front of 
the wave by a "suction" of particles from the range before the front into the 
extended range. 

The case of 5i ~ 1 (5 ~ rVn) corresponds to the smooth enough behavior 
of the deformation disturbance. Here, in the region of the front, the local 
nonequilibrium diffusion flux is nearly twice as large as that for the local equi- 
librium system. The case of 6i << 1 {S << tVb) corresponds to the sharp 
deformation variations. From figures 4 and 5 it is seen that in this case one 
can speak about the significant acceleration of the mass transfer in the local 
nonequilibrium system. The all effects noted above are by an order less in the 
local equilibrium system. 

Thus our investigation has shown that at the taken magnitudes of the parame- 
ters of the model an enhanced mass transfer takes place in local nonequilibrium 
system in comparison with the same process considered under the local equi- 
librium conditions. 

Note in conclusion that the results presented here have been obtained within 
the scope of the linear theory. The set of the hydrodynamical equations derived 
in section 2 can form the basis for the investigation of nonlinear effects. 
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Fig. 1. Normalized interstitial concentration c/Uq (normalized dimensionless diffu- 
sion flux J/ipUo) vs. dimensionless coordinate ^ at (p = 0.5, 6i = 1 and Uq > 0. 
Local nunequilibrium system - curve 1; local equilibrium system - curve 2. 

Fig. 2. Normalized interstitial concentration c/Uq (normalized dimensionless diffu- 
sion flux J/ipUo) vs. dimensionless coordinate at ip = 0.95, 6i = 1 and Uq > 0. 
Local nunequilibrium system - curve 1; local equilibrium system - curve 2. 

Fig. 3. Normalized interstitial concentration c/Uq (normalized dimensionless diffu- 
sion flux J/ipUo) vs. dimensionless coordinate ^ at ip = 1.05, di = 1 and Uq > 0. 
Local nunequilibrium system - curve 1; local equilibrium system - curve 2. 

Fig. 4. Normalized dimensionless diffusion flux J/Uq (normalized interstitial concen- 
tration ipc/Uo) vs. dimensionless coordinate ^ at <^ = 0.95, di = l/|a| and Uq > 0. 
Local nonequilibrium system - curve 1; local equilibrium - curve 2. 

Fig. 5. Normalized dimensionless diffusion flux J/Uq (normalized interstitial concen- 
tration ipc/Uo) vs. dimensionless coordinate ^ at <^ = L05, Si = l/|a| and Uq > 0. 
Local nonequilibrium system - curve 1; local equilibrium - curve 2. 

Fig. 6. Normalized dimensionless diffusion flux J/Uq (normalized interstitial concen- 
tration ifc/Uo) vs. dimensionless coordinate ^ for the case of local nonequilibrium 
system at (5i = l/|a|, L^o > and (p = 0.95 (curve 1), ip = L05 (curve 2). 
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